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Abstract. We present work in a model used to describe semi-inclusive deep inelastic scattering 
CN off the deuteron. The model uses the virtual nucleon approximation to describe the interaction 

of the photon with the bound neutron and the generalized eikonal approximation is applied to 
calculate the final-state interaction diagram. Comparison with data taken at Jefferson Lab shows 
good agreement in the covered range of kinematics and points at a largely suppressed off-shell 
\^ rescattering amplitude. The W and dependences of the total cross section and slope factor of 

the interaction of DIS products, X, off the spectator nucleon are extracted. Starting from the JLab 
"■"^ data and our model calculations, we outline and apply an extrapolation method to obtain the neutron 

structure function F2N at high Bjorken x. 
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^ INTRODUCTION 

> 

00 A reaction that can be used to study the influence of QCD dynamics at nucleonic length 

^ scales is semi-inclusive deep inelastic scattering off the deuteron (d{e^ e'ps)X), where a 

^ spectator proton is detected in the final state. At very low spectator momenta the neutron 

^ is nearly on-shell and the reaction can be used to extract information about the "free'' 

O neutron structure function in a way that minimizes the nuclear binding effects inherent 

^ in scattering off a deuteron target. At higher spectator momenta the reaction can be 

!^ used to study the modifications of nucleon properties and the role quark degrees of 

. ^ freedom in high density configurations of the deuteron. In kinematics wich favour final- 

^ state interactions (FSI), the space-time evolution of hadronization can be examined. Two 

c3 recent Jefferson Lab Hall B experiments have studied the reaction: one at high [1], and 

the other at low spectator momenta [2]. 

In order to provide meaningful interpretations of the measured data, theoretical mod- 
els that quantify the importance of the FSI in the reaction are needed. The major dif- 
ficulty in doing this is that one lacks detailed information about the composition and 
space-time evolution of the hadronic system produced in the deep inelastic scatter- 
ing and how this changes as a function of Bjorken x and Q^. Several theoretical ap- 
proaches that provide models for the d{e^e'ps)X reaction can be found in the literature 
[3, 4, 5, 6, 7, 8, 9, 10]. In the following, we will present work in a model based on the 
general properties of soft rescattering [11]. Results are compared to data taken in the 
JLab Deeps experiment and we use these results and data to describe a method to extract 
the neutron structure at high x. 



FORMALISM 



The interaction of the virtual photon with the bound nucleon is treated with the virtual 
nucleon approximation (VNA) [4, 5, 12]. In the VNA, the spectator nucleon is taken on 
the mass shell, while the virtual photon interacts with an off-shell nucleon. The VNA 
is based on the following main assumptions: (i) only the pn component of the deuteron 
wave function is considered in the reaction, (ii) the negative energy projection of the 
virtual nucleon propagator gives negligible contribution to the scattering amplitude, and 
(iii) interactions of the virtual photon with exchanged mesons is neglected. Assumptions 
(i) and (ii) can be satisfied when the momentum of the spectator proton is limited to 
Ps < 700 MeV/c [13], while (iii) is satisfied at large (> 1 GeV^) [14, 15]. 

The nuclear wave function in the VNA is normalized to account for the baryon number 
conservation [16, 17, 18, 19]: 

J a\^^D{p)\^d'p=h (1) 

where a = 2 Md' ^^S^^ ^^^^ momentum fraction of the deuteron carried 

by the bound nucleon normalized in such a way that the half of the deuteron momen- 
tum fraction corresponds to a = 1. Because of the virtuality of interacting nucleon it is 
impossible to satisfy the momentum sum rule at the same time, which can be qualita- 
tively interpreted as part of the deuteron momentum fraction being distributed to non- 
nucleonic degrees of freedom. 

In the model, two Feynman diagrams are taken into account: the plane- wave diagram 
where the on-shell proton is a pure spectator and the FSI diagram where a soft rescatter- 
ing occurs between the produced X and the spectator. To describe the FSI of the produced 
X with the spectator nucleon, the generalized eikonal approximation is applied. At the 
energies under consideration in the JLab experiments, we can assume the rescattering to 
be diffractive and occurring over small angles. The following eikonal form is adopted 
for thtX —N rescattering amplitude: 

fx'N^XN = G{Q^,XBj){i + e{Q^,XBj))e^^\ (2) 

where CJ, £ and B are the effective total cross section, real part and the slope factor of 
the diffractive X'N XN scattering amplitude. In the derivation of the FSI amplitude, 
a factorized approach is used, whereby the interaction of the virtual photon with the off- 
shell nucleon is taken out of the integration over the intermediate spectator momentum. 
Such an approximation is also applied in quasi-elastic scattering [20, 13], being valid up 
to spectator momenta around 400 MeV/c. 

We define the lab frame four-momenta of the involved particles as pu = (M/),0) for 
the deuteron, q= {y^q) for the virtual photon (with the z-axis chosen along q), ps = 

{Es = ^pj + mj,ps) for the spectator proton and p^ = {Ex.px) = {v+Mo-Es^q-ps) 
the center of mass momentum of the undetected produced hadronic system X. The initial 
momentum of the bound off-shell neutron is defined as pi = po — Ps- After applying 
all the abovementioned approximations and averaging over (the angle between the 



electron and ^Ps plane), we can write for the differential cross section [11]: 
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aEM is the fine-structure constant, —Q^ = v^ — q^ is the four-momentum transfer, m„ 



the neutron mass, y = C)Ci = aq = x = 2^ is the Bjorken x for a moving 
nucleon, Oe is the electron scattering angle, /? = ^ 0. 18 is the ratio of the longitudinal 
to transverse cross section for scattering off the nucleon, and F2N is the effective nucleon 
structure function, which is defined at x and in principle could be modified due to the 
nuclear binding (see e.g. Ref. [4]). We use the SLAC parametrization of Ref. [21] for 
Fin- 

In Eq. (3), the distorted spectral function S^{ps) contains the contributions from the 
plane-wave and FSI diagram and takes the following form: 
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where (^^^) is the undistorted (distorted) deuteron wave function with spin pro- 
jection sd, ^{s^mx') = \/{s— {rUn — mx'Y){s — {nin + mx'Y) (with s the Mandelstam 
variable of the rescattering process), rux' is the invariant mass of the produced hadronic 
state before the rescattering and p^/^ = ps^z ~ with 
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where 7 = + 2m„v — is the produced DIS mass off the stationary nucleon. The 
two regimes in Eq. (5) originate from the condition > m^,, which can be inferred 
from the approximate conservation law for the " component in high energy small 
angle scatterings [11]. 

The on-shell rescattering amplitude takes the form of Eq. (2). For the off-shell am- 
plitude /^f^ there is no clear prescription, but following our main goal of studying 
the semi-inclusive DIS based only on basic properties of the high-energy scattering we 
identify two extreme cases for off-shell part of the rescattering amplitude, one when it is 
taken to be zero (no off-shell FSI) and the other in which off-shell amplitude is assumed 



to be equal to the on- shell amplitude /^^Vxa^ referred as maximal off-shell FSI. A third 
approach is to parametrize the amplitude as 

/•off _ ^on -^{x,Q^)t /z^x 

with \i an extra parameter that can be fitted, providing a measure of the suppression of 
the off-shell amplitude. This approach is referred to fitted off-shell FSI. 

RESULTS 

Fig. 1 compares our model calculations with a selection of data taken in the Deeps 
experiment (for a definition of the spectral function P{ps, cos Gs) see Ref. [1]). For each 
Q^^W, the total cross section a and slope parameter B in Eq. (3) were taken as free 
parameters and fitted to the data for all ps values. In the case of the fitted off- shell FSI, jU 
was taken as an extra free parameter. As can be observed in the figures, the calculations 
including FSI do a fair job of describing the data over the kinematic range of the 
experiment. The differences between the three different off- shell descriptions becomes 
smaller with higher ps. At ps = 300 MeV the plane- wave and FSI amplitudes are of equal 
magnitude, making the final result very sensitive to the different descriptions, especially 
noticeable in the backward region. At this ps value, there is also an oscillating structure 
in the data which disappears with higher W, but is still present in the model results. 

At the higher spectator momenta, the calculations including FSI continously grow for 
angles in the forward direction, clearly different from the fairly flat behavior of the plane- 
wave calculations. This high FSI contribution in the forward region is a consequence 
of the structure of the phase factor in Eq. (5). At the highest ps, the calculations 
systematically underestimate the data, hinting at a breakdown of the factorization used in 
the derivation of Eq. (3). Over the whole kinematic range of the data, the fitted off-shell 
FSI calculations are more in agreement with the no off-shell than the maximal off-shell 
ones, pointing in the direction of a largely suppressed off-shell rescattering amplitude. 

The values of the fitted total cross section a and slope parameter B used in Fig. 1 are 
shown in Figs. 2 and 3 for the no off-shell FSI and maximal off-shell FSI calculations. 
After an initial peak at W = 1.2 GeV, corresponding with the creation of a A, the value 
of G drops to around 25 mb and rises with increasing W. This agrees with the picture 
of the increased creation of hadronic constituents as W increases. With increasing Q^, 
the value of the (XN) cross section parameter also becomes consistently smaller in this 
region, indicating reduced final-state interactions. This could point at the onset of a color 
transparency signal, but more data points at other and W values are needed to make 
more substantial claims. The value of the slope parameter B is also largely correlated 
with the 2^, W-dependence of the a parameter. 

NEUTRON STRUCTURE AT HIGH x 

We now use the Deeps data and our model calculations to outline a method to extract 
the neutron structure function F2 at high x from kinematics at low x. It is based on 
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FIGURE 1. (Color online) Comparison between the Deeps data [1] and model calculations. The dashed 
black curve is a plane-wave calculation, the other include final-state interactions. The effective total cross 
section and slope parameter in the final-state interaction amplitude are fitted parameters for each Q^^W, 
the real part is fixed at e = —0.5. The dot-dashed green curve only considers on-shell rescattering, the 
dotted blue curve has an off- shell rescattering amplitude equal to the on-shell one and the full red curve 
uses the off-shell parameterization of Eq. (6) 




FIGURE 2. (Color online) The fitted values of effective cross section a and slope factor B for the no 
off-shell FSI calculations used in Fig. 1 as a function of the invariant mass W. Full blue curve is for 
= 1.8 GeV^, the dashed green curve for = 2.8 GeV^. 
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FIGURE 3. (Color online) The fitted values of effective cross section a and slope factor B for the 
maximal off-shell FSI calculations used in Fig. 1 as a function of the invariant mass W. Full blue curve is 
for 2^ = 1.8 GeV^, the dashed green curve for = 2.8 GeV^. 



an analytical continuation of the amplitude to the unphysical limit — m^ — ^ 0, 
which corresponds to an on-shell neutron. Due to the low binding energy of the deuteron 
this singularity of the amplitude is very close to the physical region. It is analogous to 
the Chow-Low procedure applied in the extraction of 7i7i and NN cross sections from 
p{n^p)X and d{n^n)pn reactions [22]. It can be shown that the plane-wave part of the 
distorted spectral function in Eq. (4) has a double pole in t\ while the FSI part hasn't 
[5]. This allows us to extract the F2 structure function as 
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where j^"^ are the experimentally measured deuteron structure functions. The rhs of 
Eq. (7) has a quadratic dependence on t' . 

The approach we now take is to take a trajectory through the Deeps data at a fixed 
angle that starts at high W^ps and goes to low W^ps values. This translates into a 
descending series in ^ and an ascending series in i ^ x, as is shown in Fig. 4. 
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FIGURE 4. (Color online) Example of a trajectory int^x used in the extrapolation procedure for the 
neutron structure function F2. The z-axis shows the ratio of Eq. (7) to the parametrization of Ref. [21]. 
The blue (green) curve shows the plane- wave (ESI) model calculations along the trajectory, blue dots are 
JLab Deeps data [1]. 



In this manner, we can use measurements at low values of i - where F2N is well known 
- to extrapolate to regions with high x. 
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FIGURE 5. (Color online) Example of a extrapolation for F2N at x = 0.7. The blue (green) curve shows 
the plane- wave (ESI) model calculations along the trajectory, black dots are JLab Deeps data [1]. The 
dashed curves are a quadratic fit to the data points, the intercept at = yields the F2n{x = 0.7) value. 

Fig. 5 shows an example of a quadratic fit to the Deeps data, extrapolated to = to 




obtain the neutron structure function F2N at Bjorken x = 0.7. As the available data sits 
on one arm of the quadratic curve, it's not straightforward to do the extrapolation. Here, 
we imposed an extra constraint on the quadratic fit, that the t' value of its minimum 
has to coincide with the minimum of the FSI curve along the trajectory. The position 
of this minimum is related to the size of the FSI term in the distorted spectral function 
which the model does a good of describing for the kinematics along this trajectory. The 
three different fit curves were obtained by taking the data point and error bar values 
respectively for the lowest t' data point. With the current data set, the method doesn't 
provide a robust prediction for the high x neutron structure function, as trajectories with 
different spectator angles yield a range of values for F2N' Ideally, more data at lower ps 
values would provide better constraints and yield better predictions. 
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